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ABSTRACT 

^0 ■ The properties of accretion disks around stars and brown dwarfs in the a Ori cluster 

(-h ' (age 3 Myr) are studied based on near-infrared time series photometry supported by 

pi . mid-infrared spectral energy distributions. We monitor ~ 30 young low-mass sources 

over 8 nights in the J- and K-band using the duPont telescope at Las Campanas. 
We find three objects showing variability with J-band amplitudes 0.5 mag; five 
additional objects exhibit low-level variations. All three highly variable sources have 
been previously identified as highly variable; thus we establish the long-term nature of 
their flux changes. The lightcurves contain periodic components with timescales of ~ 
0.5 — 8 days, but have additional irregular variations superimposed - the characteristic 
' behaviour for classical T Tauri stars. Based on the colour variability, we conclude that 

hot spots are the dominant cause of the variations in two objects (#19 and #33), 
including one likely brown dwarf, with spot temperatures in the range of 6000-7000 K. 
For the third one (#2), a brown dwarf or very low mass star, inhomogenities at the 
. inner edge of the disk are the likely origin of the variability. Based on mid-infrared 

• ' data from Spitzer, we confirm that the three highly variable sources are surrounded 

by circum-(sub)-stellar disks. They show typical SEDs for T Tauri-like objects. Using 
SED models we infer an enhanced scaleheight in the disk for the object #2, which 
favours the detection of disk inhomogenities in lightcurves and is thus consistent with 
the information from variability. In the a Ori cluster, about every fifth accreting low- 
mass object shows persistent high-level photometric variability. We demonstrate that 
estimates for fundamental parameters in such objects can be significantly improved 
by determining the extent and origin of the variations. 
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1 INTRODUCTION 

The study of circumstellar disks has been a centerpiece 
of star formation research over the past decades. The 
relevance of disks is twofold: On one side, they consti- 
tute the remnant of the cold gas and dust reservoir from 
which the central star has originally formed, as well as the 
reservoir of mass and angular momentum for the newly 
formed star. Therefore, disk properties can be used to 
probe models of cloud collapse, fragme ntation, accretion , 
and angular momentum regulation (e. g., iKlein et al.ll2007l ; 
iBonnell et all 120071 ; iHerbst et ai1l2007l ). On the other side, 
disks provide the raw material for the buildup of plane- 
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tary systems and can thus be used to trace the onset of 
planet formation processes, e.g. the growth from ISM-like 
to mm-siz ed grains and the evolution of inner holes and 
gaps (e.g..lDullemond et' al. 2007; Paardekooper fc Mellemal 
12004 iRice et al.ll2006l ; lAlexander fc Armitage!l2007l ). 

The generally accepted ideal way of investigating the 
properties of circumstellar disks is high-resolution imaging 
in combination with an analysis of the full spectral en- 
ergy distribution (SED), usu ally aided by radiati ve transfer 
modeling (see the review by I Watson et al.ll2007fl . This ap- 
proach, however, is limited to the relatively few targets for 
which the disk can be resolved, and therefore works best 
for physically large disks in nearby star forming regions. 
A third and often neglected way to as sess disk properties 
is multi-filter variability analysis (e.g. iFernandez fc Eiroal 
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19961; ISkrutskie et~aH Il996l ; iLiu et all Il996l ; IWood et all 



20001 '). The presence of disks and accretion affects the fluxes 



over a broad range of wavelengths, from the UV to the mid- 
infrared. The particular advantage of multi-filter monitor- 
ing is that it allows us to probe dynamic processes in the 
disk, e.g. gas accretion and magnetic interaction between 
star and disk, as well as deviations from the commonly used 
2D thin-disk approximation, e.g. clumps, warping or spi- 
ral arms in the disk. In exceptional cases, s uch as the T 
Taur i stars KH15D (e.g . iHerbst et ail 120021 ') and AA Tau 
(e.g. iBouvier et alj [20031 ). variability analysis has led to an 
in depth picture of the inner disk region. 

Here we present a combined analysis of near- 
infrared/optical lightcurves and Spitzer photometry for 
a sample of variable very low mass (VLM) sources in 
the young open cluster around a Orionis (hereafter simply 
a Ori) . This particular cluster has a number of benefits for 
studying disk properties and evolution. In recent surveys, 
a Ori has proven to harbour a rich population of young 
stars, brown dwarfs, and isolated objects with planetary 
masses down to ~ 0. 008 Mq (e.g. 
200d : iBeiar et al.ll200ll : iKenvon et a! 
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2005; Caballcro et al 



20071. and references therein). With an age of about 3 Myr 
( Zapatero Osorio et al.ll2002l : ISherrv et alJl2008h , the cluster 



is intermediate in the evolutionary sequence between very 
young star forming regions like Chamaeleon, ONC, Taurus 
and the slightly older Sco-Cen, TWA, and /3Pic associations. 
About one third of the lo w-mass objects in this c luster are 
known to harbour a disk (|Hernandez et al"1l2007f ). and the 
disk fraction does not decline down to the lowest masses 
( Scholz & Javawardhana 2008 ). Finally, a Ori exhibits low 
interstellar extinction I Sherrv et al.|[200a ). which greatly fa- 
cilitates the interpretation of photometric data. 

This paper is structured as follows: The following Sect. 
[2] is focused on the observations and data reduction proce- 
dure. In Sect. 2] we probe for generic and periodic variability 
and look for variations on timescales ranging from hours to 
years. The characteristics of the variability are related to 
various physical scenarios in Sect. [5] According to our anal- 
ysis, large-scale photometric variability in VLM objects in 
a Ori is consistent with being caused by accretion hot spots, 
extinction changes, and variable disk emission. We support 
these findings with SED analysis for highly variably objects 
based on Spitzer photometry (Sect. [6)). The results from this 
study are discussed in Sect. [7] 



2 OBSERVATIONS AND DATA REDUCTION 

The main dataset for this paper was obtained in a near- 
infrared imaging campaign in November 2005. Using the 
WIRC camera on the 2.5-m DuPont telescope at the Las 
Campanas Observatory, we monitored a field in the a Ori 
cluster over 8 consecutive nights (13-20 November). In this 
run, we observed alternately in the near-infrared bands J 
and K, centered at wavelengths of 1.25 and 2.2 /im. The 
camera comprises four 1024 x 1024 chips, each covering a 
field of 200" x 200" . The chips are put together in a 2 x 2 
array with 175" gaps. Thus, with 4 pointings the camera 
gives a contiguous 13' x 13' field of view. 

Each time series image is produced by mosaicing 4 
pointings, see below. Except for the first night, we split the 



observing time per pointing in two single exposures, to in- 
crease the dynamic range. The exposure time per frame was 
20-40 s, chosen to avoid saturation of too many target stars 
in the field. Typically this results in 7-10 epochs per filter 
per night. Some nights were affected by cloud coverage and 
bright moon. Over the full observing run, the seeing was 
stable between 0'.'5 and / .'8. In night 4, one quarter of one 
chip became inoperational and remained dead for all subse- 
quent nights, causing gaps in the spatial coverage. In total, 
we obtained images for 65 epochs in J- and 66 epochs in 
K-band. 



The reduction is carried out using a pipeline based 
on routines in IRAF and Starlink, described in detail in 
iGilbank et al.l (|2003l ). Flatfielding and sky subtraction are 
carried out on a chip-by chip basis. In the first step, the 
raw images are multiplied with a bad pixel mask. A pre- 
liminary skyflat is produced by combining all images taken 
within a ~ 15min timespan and scaling to a median of 1.0. 
Dividing the raw images with this Hatfield yields the 'first 
pass' reduced images. Objects are detected and an object 
mask is derived. After applying this mask to the raw im- 
ages and thus eliminating the objects, an improved skyflat 
is produced using the same procedure as in the first pass. 
Flatfielding with this new skyflat finishes the 'second pass' 
reduction. 



The flatfielded images for each epoch (4 chips, 4 po- 
sitions) are coadded after matching stars in the overlap- 
ping areas of the four mosaic positions. The pixel values are 
scaled to a uniform exposure time. Finally, a WCS coordi- 
nate system is fit for the images using the wcstools package. 
The procedure yields one mosaic per filter per epoch. Note 
that the epoch is not uniform in a mosaic, because it takes 
~ lOmin to cover the four positions in one filter. Thus, we 
are not sensitive to variability on very short timescales. 

In addition, our analysis partially uses data from a sec- 
ond observing run. With IMACS at the 6.5 m Baade tele- 
scope at Las Campanas Observatory, a large field in the 
a Ori cluster was monitored over three consecutive nights 
on February 1-3 2006, i.e. 2.5 months after the near-infrared 
campaign. IMACS was used in short camera mode, with a 
total field of view of 27' x 27' . The camera has 8 CCDs in 
a 2 x 4 array and a scale of 0.2" per pixel. To block the 
light from the bright stars in the cluster, particularly from 
a Ori itself, we used a multi-object mask with 20" x 20" 
holes centered on the positions of very low mass members 
of a Ori and comparison stars. In total, we obtained 27, 24, 
51 images in nights 1-3, respectively, with exposures times 
of 300 sec and typical seeing of 0.8-1.0". 

Standard IRAF routines were used for the data reduc- 
tion of the IMACS images, including bias subtraction and 
flatfielding with domeflat exposures. We use the IMACS 
dataset to provide a c omplementary opt i cal lig htcurve for 
the object #33 from IScholz fc Eisloffell (|2004f ). which is 
found to be highly variable. Other known a Ori members 
with high variability are not covered in the IMACS run. The 
full results of this campaign will be published elsewhere. 
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Figure 1. RMS vs. relative magnitude for the J-band (upper 
panel) and K-band (lower panel) data from the WIRC run: solid 
lines show the binned median rms, crosses mark the VLM candi- 
date members from SE04 plusses the candidate members from 
ISherrv et al Highly variable objects are marked with 

squares (#2, #19, #33, and one additional source discussed in 
detail in the text). 



3 PHOTOMETRY AND RELATIVE 
CALIBRATION 

Aperture photometry for the WIRC images is performed us- 
ing customized IRAF routines developed in the framework 
of pre vious monitoring campaigns (see IScholz fc Eisloffel 
(2004) for more information). As input for the photometry, 
we used an object catalogue generated by running SExtrac- 
tor (|Bertin fc Arnouta 1 199a ) on a reference image. We de- 
termined pixel offsets with respect to the reference frame for 
the whole set of images and used them to produce individ- 
ual object catalogues for each time series image. We choose 
a constant aperture radius of 8pix (l'/6) for the photometry, 
which is approximately twice the radius of the average PSF. 
For sky substraction we use an annulus with 12 pix (2'.'4) 
inner and 17 pix (3'.'4) outer radius. 

The time series photometry is calibrated using a stan- 
dard lightcurve generated from the data for a set of non- 
variable reference stars. This standard lightcurve reproduces 
the effects of variable seeing and atmospheric extinction. To 
select the reference stars, we first consider only objects with 
photometry errors < 0.1 mag in all our time series images. 
To remove the variable stars from this initial set, we compare 



the lightcurve of each object to the average lightcurve of all 
other objects in the set. For each object, we subtract the av- 
erage lightcurve and calculate the rms (root mean square) of 
these deviations. Objects with large rms are considered vari- 
able and removed from the set until only the non- variables 
remain. For more details on this procedure see SE04. The 
final standard lightcurve is calculated as the average time 
series of ~ 50 reference stars in J and ~ 30 in K. We then 
subtract the standard lightcurve from that of every object 
to complete the relative calibration. 

In Fig. [T]we demonstrate the quality of the time series 
photometry by plotting the rms vs. relative magnitude. For 
this plot, we excluded 3<r outliers and all objects with fewer 
than 60 datapoints after excluding the outliers. The solid 
lines mark the median rms and thus give an estimate for 
the photometric noise in the lightcurves. As can be seen 
in the figure, we achieve an optimum accuracy of 1-2% in 
both bands. For the non- variable reference stars used in the 
relative calibration, we obtain an average rms of 0.017 mag 
in J- and 0.014 mag in K-band. Over four magnitudes, the 
noise is < 0.03 mag in J and < 0.04 mag in K. Note that 
the substellar mass limit for the a Ori cluster corresponds 
to a relative magnitude of ~ 0.0 in J and ~ 0.3 in K, i.e. 
the photometric noise is in the 2% range down to the brown 
dwarf regime. 

Photometry and calibration for the IMACS run followed 
the same recipes. Aperture photometry was obtained for 
our target #33 and comparison stars on the same chip, us- 
ing a fixed aperture of 12 pixels. For relative calibration, 
we used the average lightcurve of six comparison stars. Af- 
ter subtracting this reference time series, the comparison 
lightcurves show a rms of 0.015 mag and no clear signs of 
variability. The source #33, on the other hand, clearly ex- 
hibits large-scale variations, which will be discussed in detail 
below. 



4 VARIABILITY IN a ORI MEMBERS 

The pointing of the WIRC observations was chosen in a way 
to maximize the coverage for the sample o f clust er mem- 
ber candidates selected bv lScholz fc Eisloffell l|2004 SE04 in 
the following) based on photometry in five bands (RIJHK). 
The timeseries cover 30 out of 135 objects listed in SE04, 
from which 21 are unsaturated in all J- and K-band im- 
ages. These sources, marked with crosses in Fig. [TJ are the 
primary sample for the variability analysis. They cover a 
magnitude range of 12-17 in J (11-15 in K), corresponding 
to approximate masses of 0.02 to 0.6 M©, assuming a dis- 
tance of 420 pc (|Sherrv et al.| [2~008) and negligible reddening. 
In addition, our field co vers 12 objects from the photomet- 
ric sample published by ISherrv et al.l ([2004) , 9 of them are 
unsaturated in both bands, from which 5 coincide with SE 



4.1 Probing for generic variability 

In a first step, we probe for generic variability by compar- 
ing the rms values for the candidates with the median rms 
derived from field stars. If the rms in a particular lightcurve 
is significantly higher than the median rms at the same rel- 
ative magnitude, the object is considered variable. For the 



4 Scholz et al. 



comparison we use the F-test, with F = 0" 2 /°mcd (c: rms 
for the candidate; o" mc d: median rms). An F- value > 2 cor- 
responds to a significance level of 99.9%. According to this 
test, seven objects from the SE04 sample are significantly 
variable in both bands (#2, #8, #17, #19, #23, #29, #33). 
The object with the running number #217 from the Sherry 
et al. sample, not contained in the SE04 catalogue, passes 
the same test. 

Three of these objects stand out based on their 
lightcurve characteristics: the SE04 candidates #2, #19, 
and #33. They are marked in Fig. \T\ with squares; their 
lightcurves are plotted in Fig. [2j see Table [1] for a sum- 
mary of their lightcurve properties. All four show photo- 
metric amplitudes > 0.5 mag in the J-band and > 0.2 mag 
in the K-band. In Fig. [1] their datapoints are clearly offset 
from the median rms in both bands. The lightcurves exhibit 
variability on a range of timescales, including continuous 
variations over days and rapid changes over hours. All three 
have been identified in SE04 as objects with large ampli- 
tude variability in the I-band. Objects no. #2 and #33 in 
particular have been observed to be highly variable in two 
observing runs in January 2001 and December 2001. With 
our new data from November 2005 we extend the time base- 
line. The large-amplitude variability clearly is of long-term 
nature, sustained over almost 5 years. 

The complementary Fband lightcurve for #33, obtained 
in Februar 2006, shows the same type of variations with a 
total amplitude of 0.42 mag and a rms of 0.12 mag (see Fig. 

clearly higher than in the references stars in the same 
field. 

A fourth object in our field of view shows a lightcurve 
comparable with the three highly variable objects discussed 
above, marked as well with a square in Fig. [2] This source is 
seen as close companion (separation 2.7") to the brighter ob- 
ject J053858 . 32-021609, which is identified as young cluster 
member by ISherrv et al. (2004) based on optical colours. 
The bright source J0538-0216 is saturated in most of our 
images. It was flagged as v ariable in the CIDA variability 
survey (|Bricefio et al.ll2005h . which covers ti mescales from 
days t o 5 years. The I-band value provided bv lSherrv et alj 
|2004h from 1996-98 is 0.8 mag fainter than the one measured 
by SE04 in 2001, another indication for variability. Since 
our photometry for the companion is certainly affected by 
light from the primary, it is likely that variability from the 
primary causes the large variations seen in the companion 
lightcurve. Lacking an information on cluster membership 
for the companion, we do not discuss this object any fur- 
ther. 

Based on the variability characteristics, near-infrared 
colours, and (for #2 and #33) strong emission lines in op- 
tical spectra, SE04 conclude that ongoing accretion is the 
most likely reason for the strong photometric variations. In 
Sect. [5] we will further explore this scenario based on our 
near- infrared lightcurves. 

The remaining five variable objects show an rms signifi- 
cantly increased in comparison with the median rms, in both 
bands. In these cases, the total photometric variations are 
< 0.2 mag; the rms < 0.05 mag. Such low-level variability is 
seen in a number of objects in SE04, including #8 and #23 
(but not #17 and #29), and is attributed to the presence of 
cool spots, co-rotating with the objects. The rms for all five 



Table 1. Object and lightcurve properties for the four highly 
variable objects. The average J- and K-band magnitudes in the 
first two rows are derived following the calibration given in Sect. 
14.31 Masses and effective temperatures are taken from SE04. For 
#19, the effective temperature is estimated from the photometry. 
Given the significant variability in fluxes and colours, these values 
are subject to large uncertainties (see the note in Sect. l4~3ll . 





#2 


#19 


#33 


mass (M©) 


0.07 


0.65 


0.10 


T eff (K) 


2800 


~ 4000 


3000 


< J > 


14.87 


13.15 


14.75 


<K> 


13.72 


11.42 


13.55 


J ampl. (min-max) 


0.50 


0.55 


0.91 


K ampl. (min-max) 


0.22 


0.33 


0.68 


J rms 


0.11 


0.15 


0.24 
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0.10 
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Figure 3. I-band lightcurve for object #33 obtained with IMACS 
in the three nights Feb 1-3 2006. 



objects are similar in J- and K-band (withi n ±0.01 mag), as 
expected for cool spots (|Scholz et al.ll2005r ). 

In combination with literature results, there is now a 
sample of 4 a Ori objects with masses close to or below 
the substellar limit, showing high-level variability and spec- 
troscopic signature of accretion: two in this study (#2 and 
#33) , one more in SE04 (the sourc e #43), and the objec t 
SOriJ053825.4-024241 discussed bv lCaballero et all (|2006h . 
All four have a similar type of lightcurve, consistent with the 
typical variability signature of accretin g young stars (e.g. 
iHerbst et all 1 19941 : iBouvier et"aH Il995l ). making them the 
very low mass representatives of classical T Tauri stars. 



4.2 Period search 

If photometric variability is due to spots on the surface of 
the star, the flux is expected to be modulated periodically by 
the rotation of the object (see Sect.JSJ. Therefore, a period 
search procedure was carried out for the eight objects iden- 
tified as being variable in Sect. 14. ll One main intention here 
is to verify the periods published in SE04, in particular for 
the highly variable and accreting targets. For these objects, 
superimposed irregular variations, e.g. due to a variable ac- 
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Figure 2. Lightcurves in J-band (upper row) and K-band (lower row) for objects #2 (1st column), #19 (2nd column), and #33 (3rd 
column). The time axis is not continuous; for viewing purposes we subtract 0.5 d from the 2nd night observing times, 1.5 d from the 3rd 
night, and so on. One tickmark on the x-axis thus corresponds to two nights. Nights are separated by horizontal dashed lines. 



cretion rate, may hamper the detection of a photometric 
period, particularly if the irregularities occur on timescales 
shorter than the rotation period. In these cases, repeated 
monitoring is required to disentangle the various contribu- 
tions to the flux changes and to establish the most likely 
rotation period. 

The period search is based on the J-band data: As 
shown in Sect. [S] spots will cause more variability in J- 
than in K-band. Our procedure combines three independent 
tests, which are introduced in the following; for more d etails 
on these algorithms see SE04. IScholz fc Eisloffel (|2005l ). and 
the references therein. The lightcurves are prepared for the 
period search with a simple cr-clipping algorithm to exclude 
outliers. 

a) W e start the p eriod search with a Scargle peri- 
odogram (IScarglelll98"a ). claimed to retain the exponential 
distribution of peak heights for unevenly sampled datapoints 
(i.e. the probabilit y to find a peak is P(z) = exp (z) with 
z the peak height). Iflorne fc Baliunasl (|l986l ) give a empiri- 
cally found equation to estimate the False Alarm Probabil- 
ity (FAP) from the peak height in the Scargle periodogram, 
which has be en used previously to set thresholds for period 
searches (e.g. lLamm et al.ll2005D . 

b) For a robust estimate of the FAP for a given pe- 
riod, we fit the lightcurve with a sine function using a least- 
square procedure, with with amplitude and zeropoint as free 
parameters. The variance of the original lightcurve and the 
variance after subtraction of the sinewave are compared sta- 
tistically using the F-test. If the period is not significant, its 
subtraction will not alter the noise in the lightcurve in a sig- 
nificant way (jScholz fc Eisloffelll2004l ). The resulting phase- 
plots, shown in Fig. [4] for the best periods, also allow us to 
check the periods visually. 

c) The CLEAN algorithm bv lRoberts et all (|l987l ) de- 
convolves the original periodogram and the window function 
and thus 'cleans' the spectrum from sidelobes and aliases. 



The procedure can therefore be used to distinguish between 
spurious features and artefacts in the periodogram. We use 
CLEAN to check if the frequencies found in the Scargle pe- 
riodogram are artefacts. The peak heights in the CLEANed 
periodogram are not easily transformed into probabilities, 
therefore we do not attempt to use this periodogram to es- 
timate FAPs. 

In Fig. |3] we provide the datapoints as function of 
phase for the best-fitting periods. Additionally, we show the 
CLEANed periodograms for the highly variable objects in 
Fig. [S] In the following, we discuss the results for each ob- 
ject. 

For object #2, the Scargle and CLEAN periodograms 
indicate three possible periods of 13.7 h, 25.9 h, and 180 h; 
the first one agrees reasonably well with the period of 14.7 h 
given by SE04. In all three cases, however, the FAP from 
the F-test is relatively high (~ 2%). The rotation period 
may be either the ~ 14 h period seen in all three available 
lightcurves or the period of 7.5 d seen in the current dataset. 
While this long periodicity is not seen in the SE04 dataset, 
it may have easily been masked by irregular variations on 
shorter timescales. The two longer plausible periods of 25.9 
and 180 h are close to multiples of the 14 h period and may 
not be real. 

For object #19, our procedure indicates a significant pe- 
riod of 31.3 h. From the peak HWHM in the CLEANed pe- 
riodogram we estimate a period uncertainty of ±3h, which 
indicates that our period and the one in SE04 (39.3 ± 2.3 h) 
are not too different. Thus, a rotation period of 30-40 h is 
consistent with all presently available lightcurves. 

For object #33 our tests indicate two possible periods of 
34.3 h and 168 h. Applying the same period search procedure 
to the I-band lightcurve for this target gives a likely period 
of 44 h, but the corresponding peak in Scargle and CLEAN 
periodogram is broad, indicating an uncertainty in the range 
of ±10 h. Thus, this value is consistent with the shorter J- 
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band period. The optical lightcurve does not allow to probe 
the longer J-band period, due to its limited time coverage. 
In SE04 two possible periods are listed for this object, one 
is 44.6 h, the other one is 228±86h. The best interpretation 
of these findings from four monitoring campaigns is that the 
long-term lightcurve is showing evidence for two periods of 
35-45 h and ~ 200 h. At this point, it is difficult to establish 
which one is the rotation period. 

Thus, none of the three highly variable objects exhibits 
an unambiguous photometric period. Combining multi- 
season datasets, however, gives a good idea of persistent 
periodicities in the flux changes of these objects, allowing 
us to establish the most likely values for the rotation peri- 
ods, as shown in Fig. [4] All three objects show strongly non- 
sinusoidal lightcurves, i.e. significant residuals remain after 
subtraction of a sinewave with the detected period. There is 
clearly additional variability on timescales different from the 
assumed rotatio n period, as expected for accreting T Tauri- 
like objects ("e.g. lHerbst et al.lll994l ; lBouvier et alj(l995l ). As 
pointed out earlier, the confidence level for the given peri- 
ods is not easily estimated from one dataset alone, and can 
only be assessed based on repeated monitoring. The periods 
published here should still be considered preliminary and 
require further verification. 

For object #17 we find a period of 60.6 h with an am- 
plitude of about 0.15 mag, significant in the Scargle and 
CLEAN periodogram, as well as in the F-test. This period 
may correspond to the rotation period for this object (see 
Sect. [5}, although the phaseplot does not look entirely con- 
vincing (Fig. El). The time series for #8, found to have a 
period of 79.1 h in SE04, shows no periodic signal, except 
for beat periods at 0.1, 1.1, 2.1 d. The same result is found 
for #23 and the Sherry et al. object #217. Visual inspection 
reveals a slow trend in these three lightcurves, which may 
indicate a rotation period longer than our time coverage. 
Finally, for object #29, our lightcurve does not contain any 
significant periodic signal. 

4.3 Comparison with 2MASS 

Comparing the near-infrared magnitudes from the 2MASS 
catalogue with our lightcurves allows us to carry out an- 
other test for variability on timescales of ~ 7 years. Based 
on a sample of ~ 400 field stars, we determined the off- 
sets between our relative magnitudes and the 2MASS abso- 
lute photometry to be A J = 14.35 ± 0.13 mag and AK = 
13.09 ± 0.14 mag. Most of the uncertainty in these terms 
comes from chip-to-chip variations in the zeropoint. This cal- 
ibration is reliable for objects in the range J = 12 ... 16 mag 
and K — 11 ... 15 mag; brighter sources are affected by satu- 
ration in our data, fainter ones have poor s/n in 2MASS. We 
find that for almost all objects in the a Ori sample the aver- 
age magnitude from our lightcurves agrees with the 2MASS 
magnitude within 2a. Thus, for the majority of the a Ori 
cluster members the fluxes are stable on timescales of years. 

The only two exceptions are objects #2 and #19, al- 
ready identified as highly variable (SE04, this paper); both 
objects are significantly fainter in our dataset compared with 
2MASS. The magnitude difference in J (K) between our 
lightcurves and 2MASS is 1.30 mag (0.83) for object #19 
and 0.39 mag (0.35) for object #33. This is more evidence 
for high-level variability sustained on timescales of years, as 



already pointed out in Sect. |4~T1 These two examples illus- 
trate that single epoch photometry/spectroscopy is of lim- 
ited use when deriving parameters for accreting stars and 
brown dwarfs. Large variability in fluxes and colours has to 
be taken into account and adds substantial uncertainties. 
This issue will further be discussed in Sect. [7] 



5 ORIGIN OF THE VARIATIONS 

The most commonly discussed origins of photometric vari- 
ability in young stellar objects include a) cool spots induced 
by magnetic activity, b) hot spots formed by the accretion 
flow, c) variable circumstellar extinction, and d) variable 
disk emission. Options b-d are associated with the presence 
of an accretion disk. In the following, we will provide a more 
in depth discussion of the characteristic variability expected 
for each of these scenarios. Subsequently, the predictions will 
be compared with the properties of the lightcurves (Sect. 
153) . 



5.1 Cool spots 

Cool spots, comparable to the spots on the Sun, are one of 
the most common sources of photometric variability in low- 
mass stars at all ages. If the spot distribution is not fully 
symmetric, a periodic flux modulation is caused as the spots 
co- rotate with the object. Spot activity is produced by the 
interaction of stellar magnetic fields with the photospheric 
gas, and is thus an indicator of magnetic activity. 

Typically, young stars exhibit spot filling factors (frac- 
tion of one hemisphere covered by the spots) of 0-30%. 
The temperature contrast bet ween spots and pho tosphere 
is in the range of 10-30% (|Bouvier et al.l Il995h . There 
are indications that the spot properties change in the 
very low mass re gime, possibly to lower filling factors 
l|Scholz et al.l [2005). Spot configurations typically remain 
constant over timescales of s everal days, but undergo 
changes on timescales of weeks frTussain 2003). 

With these properties, cool spots are expected to pro- 
duce strictly periodic lightcurves in our 8-day observing win- 
dow. The amplitudes of the variability are expected to be 
< .15 mag in J-band and < 0.1 mag in K-band (see Fig. 
3 in lScholz et al.|[2005l ). The objects become redder as they 
become fainter, with a J — K colour variability < 0.05 mag. 



5.2 Hot spots 

Hot spots in young stars are thought to be caused by 
infalling gas and are thus a direct consequence of accre- 
tion. Similarly to their cool counterparts, hot spots are co- 
rotating with the objects and may thus induce periodic vari- 
ability on timescales of the star's rotation period. Changes 
in accretion rate or accretion flow geometry can produce ir- 
regular flux modulations on timescales ranging from hours 
to years. As a consequence, the optical and near-infrared 
lightcurves of accreting T Tauri stars often show a complex 
behaviour. 

The photometric amplitudes expected from hot spots 
can be estimated as follows: 
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Figure 4. Phaseplots for objects with significant photometric period. For objects #2 and #19 two possible periods are shown. 
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Figure 5. CLEANed periodograms for the three highly variable objects. The frequencies coinciding with significant periods in the Scargle 
pcriodogram as well as in the F-test (see text) are marked with dashed lines. 
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Here Ts corresponds to the spot temperature and / to the 
filling factor. For this simple model, we approximate the 
fluxes F(T) with a blackbody functi on. This approach is 
simi lar to the one used, for e xample, bv lBouvier et all (|l995l ) 
and ICarpenter et alJ (|200ll ). 

Hot spots generally produce large-scale photometric 
variability, with amplitudes declining towards longer wave- 
lengths. We calculated J- and K-band amplitudes for effec- 
tive temperatures of 3000 and 4000 K and a wide range of fill- 
ing factors (2-40%) and spot temperatures (5000-12000 K), 
but limited the choice of parameters to allow only J-band 
amplitudes between 0.4 and 1.0 mag, as observed for our 
highly variable objects (Table[TJ. In most plausible cases, the 
colour changes in J — K exceed 0.1 mag, and are thus clearly 
more pronounced than for cool spots. The relative variability 
in J — K colour, measured as the quotient Aj-k/Aj ranges 
from 0.25 to 0.45 for T e g = 3000 K and from 0.15 to 0.35 for 
T c ff = 4000 K. Larger relative colour changes would require 
very small filling factors (< 1%) to produce the observed 
J-band amplitudes. 



5.3 Extinction 

Variable extinction is caused by temporally changing inho- 
mogenities in the absorbing material along the line of sight. 
For young stars, the line of sight intercepts parts of the ac- 
cretion disk as well as the surrounding molecular cloud. For 
the cloud, we expect variations on timescales of weeks or 
longer, based on the typical velocity field of the molecu- 
lar gas l|Carpenter et al . 2001). The timescale for extinction 
variations originating from the disk is determined by the ro- 
tational velocity in the disk and can range from hours to 
years, depending on the radial distance from the star. 

Possible scenarios for inhomogenities in circumstellar 
disks include azimuthal asymmetries (e.g. spiral arms), a 
partially flared outer disk, a warped inner disk, or dust i n 
an inhomogenuous disk wind (|Tambovtseva fc GrininfeoOcf ) . 
If the features are restricted to small scaleheights, the ori- 
entation of the disk would have to be close to edge-on to 
produce extinction variations. On the other hand, a face-on 
geometry is not expected to cause flux changes due to vary- 
ing extinction. Depending on the stability of the feature, the 
lightcurve may exhibit a periodic component. 

Similarly to cool and hot spots, variable extinction will 
make the objects redder as they become fainter. The wave- 
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length dependence of extinction in the near-infrared is usu- 
ally assume d to be a powe r law A x oc \~ ri . With f3 ~ 1.7, 
as used by iMathid |l990h . we obtain Aj -k/Aj = 0.62 
Lowe r values for j3, as seen in some clouds (|Froebrich et al.l 
120051 ). would lead to a steeper extinction law and a smaller 
value for Aj-k/Aj. For example, = 1.5 or 1.3 yield 
Aj-k/Aj — 0.56 or 0.52, respectively. In any case, the slope 
of the colour changes induced by variable extinction is larger 
than the values expected for hot spots and can thus be used 
to distinguish between both options. 



5.4 Disk Emission 

The dust in circumstellar disks absorbs incoming radiation 
from the star and re-radiates it at longer wavelengths. This 
produces the infrared excess emission typically observed in 
young stellar objects at A > 1 — 2 fim. Near-infrared dust 
emission strongly depends on the temperature of the emit- 
ting material. The temperature structure in the inner disk 
can be affected by the release of energy, for example due to 
changes in the accretion rate or the inner hole size. Such 
processes can thus induce near-infrared variability. 

Disk emission becomes stronger at longer wavelengths, 
i.e. we expect more pronounced variations in the K-band 
compared with the J-band. Thus, variable disk emission is 
expected to make the star bluer as it becomes fainter, in con- 
trast to all aforementioned origins of variability. It is difficult 
to set firm limits on the amplitudes expected for this type of 
vari ation. Based on i nner disk models provided by N. Cal- 
vet. lCarpenter et al.l (|200ll ) conclude that J-band variability 
up to 1 mag and J — K colour changes of a few tenths of a 
magnitude are conceivable. Changes in disk emission may 
occur on a variety of timescales, ranging from days to years. 



5.5 Comparison with the lightcurves 

The derived lightcurve characteristics for the various sus- 
pected sources of variability allow us to constrain the na- 
ture of the variations seen in our observing campaign. For 
the three highly variable sources, #2, #19, and #33, it is 
clear that cool spots can be excluded as causes of variabil- 
ity, simply because the amplitudes are too large. In all three 
cases, the amplitude in J-band is significantly larger than 
in K-band, i.e. they are redder when they are fainter, which 
excludes disk emission as the dominant source of variability. 
To decide between the two remaining options, hot spots and 
extinction, we compare the colour variability in J — K with 
the predictions derived in Sect. I5.2l and [5731 

In Fig. [6] we provide J vs. J — K plots for all three ob- 
jects. With solid lines we show linear least-square fits to the 
data. Large symbols indicate the average values for each in- 
dividual night. Dashed/dash-dotted and dotted lines show 
for models for hot spots and extinction, respectively, plotted 
with arbitrary offsets for clarity. As explained in Sect. [2] J- 
and K-band observations were obtained alternately; there- 
fore the epochs for J-band and K-band differ typically by 
15-20 min. As a consequence, the J — K value corresponding 
to a given J-band epoch is only approximately correct. This 
effect will cause additional scatter in the J vs. J — K plots, 
particularly for objects with rapid variability. In the follow- 
ing, we discuss all three highly variable objects separately. 
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Figure 6. Colour variability in comparison with models for hot 
spots (dashed and dash-dotted lines) and extinction (dotted line). 
Models are offset for clarity. Solid lines are a linear fit to the 
datapoints. Large dots show the photometry averaged for each 
night. 



#2: The datapoints follow a straight line in the dia- 
gram, with a slope of Aj_k/Aj = 0.55 ± 0.03. If we fit the 
nightly averages, we obtain a very similar value, indicating 
that variations on timescales of hours and days are caused by 
the same process. In Fig. [6] we compare with the expected 
slope for spot temperatures of 5000 and 6000 K (dashed, 
dash-dotted line) and variable extinction (dotted line) . This 
slope is inconsistent on a 3a level with the predictions for 
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hot spots, no matter what spot parameters are chosen. The 
difference between J- and K-band amplitude is simply too 
large to be explained by hot spots. On the other hand, the 
distribution of datapoints agrees well with the expectations 
for variable extinction. The total J-band amplitude in our 
observations corresponds to a change in the extinction of 
AAv ~ 2 mag. As discussed in Sect. 14.21 the lightcurve con- 
tains periodic components with timescales of 1-8 d, which 
favours the disk as the origin of the extinction changes. As- 
suming Keplerian rotation, this translates to a distance be- 
tween the star and the source of variable extinction between 

0. 01-0.04 AU, which is in the ran ge of the typical du st sub- 
limation radius for VLM objects (|Scholz et afll2006h . Thus, 
if our interpretation is correct, we are observing a feature at 
the inner edge of the disk. 

#19: The variability characteristic for this object 
presents the most difficult case in the interpretation. Strong 
scatter is seen particularly in the left/blue part of the 
lightcurve. The datapoints for nights 3-8 agree well with a 
linear fit with Aj-k/Aj = 0.41±0.02. A similar slope is ob- 
tained by fitting the averages. This value is clearly inconsis- 
tent with variable extinction, but can roughly be reproduced 
by hot spots, which are likely the dominant source of vari- 
ability for this object. In Fig. [6]we show the slopes for spot 
temperatures of 6000 K (dashed) and 7000 K (dash-dotted 
line), for comparison. In the first two nights, however, the 
object is too faint and blue for the hot spot model. As can 
be seen in Fig. [2] the source becomes gradually brighter and 
redder over the first three nights, a typical sign for variable 
disk emission. In our lightcurves, we might see the effects of 
an increase in the accretion rate, leading a) to more heating 
in the inner disk and thus enhanced disk emission (nights 1- 
3) and b) to the formation of a hot spot close to the surface 
of the star, co-rotating with the object and thus modulating 
the flux (nights 3-8). 

#33: The J vs. J — K data are well-approximated by a 
linear fit with a slope of A./_a"/Aj = 0.36±0.02. The scatter 
seen around this fit can be fully explained by a combination 
of photometric errors and the epoch difference between J- 
and K-band data. Fitting the nightly averages yields a very 
similar value. The slope is clearly too small to be consistent 
with variable extinction, as seen in Fig. [6] Hot spots pro- 
vide the most plausible explanation for the variability seen 
for this object. The measured value for Aj_x/Aj is best 
reproduced by a spot temperature of 6000 K. Our J-band 
amplitude then constrains the filling factor to > 20%. This 
should be treated as a lower limit, because we cannot be cer- 
tain that our observing window covers the total variability 
amplitude. The datapoints from two nights (nights 2 and 
5) seem to be separated from all others, clustered around 
J = 0.1 and J — K = —0.15. In these nights, the object 
reaches its maximum in brightness and the bluest colour, 
indicating that we may be seeing a hot region projected 
against the cool surface of the star. 

Finally, we note that the J vs. J — K datapoints for 
object #17, found to have a period of 60.6 h in the J-band 
data, mostly scatter within ±0.05 mag of the mean. The 
colours tend to become redder as the object becomes fainter, 

1. e. cool spots are the most likely agent of the variability. 
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Figure 7. Spectral energy distributions for the five highly vari- 
able sources in the a Ori cluster. The crosses show the SEDs for 
#2, #19, #33, #43, discussed in this paper and by SE04; the 
square shows the SEP for SCMJ053825.4-024241, discussed by 
ICaballero et all (I2006T) ■ All fluxes are scaled to the J-band flux. 
For reference, overplotted as dots are the datapoints for all objects 
with disks in the SE04 sample, see Appendix lAl 



6 SPECTRAL ENERGY DISTRIBUTIONS 

As the origin of the high-level variability is assumed to be 
either hot spots or clumpy circumstellar material, and thus 
associated with the accretion disks, we aim to improve the 
variability analysis by investigating the near/mid- infrared 
spectral energy distributions (SED) for the highly variable 
sources in the a Ori cluster. Spitzer IRAC and MIPS pho- 
tometry was carried out for the objects #2, #19, #33, #43 
discussed here and in SE04, see Appendix[X]for information 
on the photometry. For the variable sourc e SOriJ053825.4- 
024241 identified bv lCaballero et al.l (l2006t). we obta i ned th e 
Spitzer data from the survey by Hernandez et al.l l|2007h . 
which classifies the source as Class II object. 

The SEDs are complemented by near-in frared data from 
2MAS S and I-band data from either SE04 or lCaballero et al.1 
(2006). It should be noted that the optical, near-infrared, 
and mid-infrared data are taken at different epochs. The 
SEDs are thus expected to be affected by variability. This 
particularly applies to the optical and near-infrared bands 
where the variations are strongest (at least for #2, #19, 
#33, and SOriJ053825.4-024241; not necessarily for #43). 

The fluxes for the five highly variable objects are given 
in Table [5] We show their SEDs, scaled to the J-band flux, 
in Fig. [7] For comparison, we overplot the SEDs of the total 
disk sample in SE04 (see Appendix [XJ. All five highly vari- 
able objects have SEDs typical for classical T Tauri stars. 
None of them shows an extreme or highly unusual SED in 
comparison with the total disk sample. (This is also appar- 
ent from the colour plots given in Appendix^]) Thus, the 
presence of variability is not necessarily correlated with un- 
usual features in the SED. 

The four objects from SE04 are relatively similar in the 
IRAC wavelength regime (3.6-8 /im). At 24 /im, however, 
two of them, #2 and #43, have datapoints around -1.0, 
while the other two are about 0.6 dex below that. To fur- 
ther explore the data, we use model SEDs generated from a 
Monte Carlo radiation transfer code, the same set of models 
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Table 2. Photometric data for the five hi ghly variable sources in 
the a Ori cluster. I-band data comes from lCaballero et al.l (2006) 
for SOriJ053825. 4-024241 or from SE04 for all remaining objects. 
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as applied in lScholz et al.l (|2006l ). In the following, we give a 
brief description of the main features of the code, for further 
information see the aforementioned paper. 

The models use NextGen model atmospheres (e.g. 
lHauschildt et all Il999l ) for the photospheric spectra, with 
log g = 4.0. We apply a distribution of grain sizes following a 
power law with an exponential decay for particles with sizes 
above 50/xm and a f ormal maximum grain size of 1 mm, see 
I Wood et all (l2002bl ) for a more detailed discussion of the 
dust grain model. For all models we assume that dust in re- 
gions close to the star is destroyed if temperatures rise above 
the value for dust sublimation, which sets a minimum inner 
dust radius of ~ 7i?*. Within that radius, the disks are as- 
sumed to be optically thin. In the models, the scaleheight 
of the disk increases with radius, h(r) = ho (r/R*) 13 , so that 
the degree of flaring can be varied by adjusting f3 and ho- 

Fitting SED da ta is subject to a num ber of degeneracies, 
see the discussion in lChiang et alj l|200ll ). T ypically, a range 
of mo dels is able to fit the observed SEDs (|Robitaille et al.l 
2007). Therefore, we do not claim that the models presented 
here are a unique solution for the disk structure. Instead, 
they should be seen as exemplary fits to the SED. For the ob- 
jects considered here, further complications are introduced 
by the significant variability in the optical and near-infrared 
bands. 

To facilita te the modeli n g, we fix a number of parame- 
ters. Following ISherrv et all l|2008l ). the average interstellar 
extinction is assumed to be Ay — 0.5 mag and the distance 
420 pc. Since the near- and mid-infrared SED is not sensi- 
tive to global disk parameters, we fix the total disk mass 
at 10~ 3 Mq and the outer disk radius at 100 AU. a Ori is a 
cluster without evidence for deeply embedded sources, there- 
fore we fit the SEDs without including an envelope. Addi- 
tionally, we assume that accretion rates are fairly low and 
disk heating is thus dominated by radiation from the central 
sources. With these assumptions, the main free parameters 
are the flaring power, the inclination, and the inner disk 
radius (jWood et al.ll2002al ). 

To explore possible reasons for the discrepancy in the 
24 fim flux, we match the SEDs of #2 and #33, two ob- 
jects with similar effective temperature, with models. We 
find that only one parameter, the flaring power, needs to be 
modified to account for the difference. This is illustrated in 
the two panels in Fig. [5] The best fit of the optical/near- 
infrared data is produced with T c ff = 2900 K for object #2 
and 3000 K for object #33. For the two disks we use essen- 



tially the same model parameters, with the exception of the 
flaring power. We assume ho = 0.025, j3 = 1.2 for #2 and 
h = 0.02 and (3 = 1.1 for #33. At 1 AU distance from the 
star, this results in a scaleheight of 29 i?* for object #2 and 
10 7?* for object #33. As seen in the plots, such a choice of 
parameters can nicely explain the 0.6 dex flux difference at 
24 fim in these two SEDs. The larger scaleheight leads to 
a larger surface area that is irradiated by the stellar light 
and thus an increased level of mid-infrared flux. Similarly, 
the SED of source #43 can be explained with an enhanced 
flaring power. 

This provides additional constraints on the scenario pre- 
sented in Sect. !5.5l to explain the variability. While the pho- 
tometric variations for #33 and #19 are best explained by 
hot spots, object #2 likely undergoes changes in circum- 
stellar extinction, caused by inhomogenities at the inner 
disk edge (see Sect. 15. 5p . For this to have an effect on the 
optical/near- infrared fluxes, the line of sight has to intercept 
parts of the accretion disk (Sect. [53)) . A disk with compar- 
atively high degree of flaring and thus relatively large scale- 
heights would naturally comply with this prerequisite. On 
the other hand, relatively small scaleheights favour the pos- 
sibility of having a unobstructed line of sight towards the 
hot spots in the inner accretion zone, as it is the case for 
#19 and #33. Hence, the results from the SED modeling fit 
into the interpretation given in Sect. 15.51 



7 DISCUSSION 

In Sect. |4l we have identified three low-mass objects in 
a Ori showing strong variability with amplitudes exceeding 
0.5 mag in the J-band, two of them with masses close to or 
below the substellar limit. Their lightcurves are not strictly 
periodic, instead we see evidence for multiple periods with 
additional irregular contributions. A few more likely VLM 
objects in a Ori exhibit small-scale photometric variations 
with amplitudes < 0.1 mag. 

This type of variability is in line with the phenomeno- 
log ical classificat i on of variability in T Tauri stars, as given 
bv lHerbstetaD dHH): The low-level modulations can be 
classified as 'type F variations, mostly seen in diskless stars, 
while the large-scale, partly irregular modulations resemble 
the 'type II' variability, predominantly seen in disk-bearing 
objects. Thus, the variability characteristic of very low mass 
stars and brown dwarfs is completely analogous to the one 
seen in solar-mass T Tauri stars. 

For two highly variable objects, the colour variabil- 
ity is consistent with hot spots, formed by the shockfront 
of the accretion flow. Since hot spots are co-rotating with 
the objects, they are expected to induce a periodic flux 
modulation; additionally, changes in the flow properties 
can cause irregular variations. For object #33, possibly a 
brown dwarf, the best fit spot temperature is ~ 6000 K 
with a filling factor of > 20%. While this value is simi- 
lar to the cons traints for hot spots in othe r brown dwarfs 
(6500- 8000 K, iHerczeg fc Hillenbrand! (|2008h see also iKoenl 
it differs considerably from assumptions made in 
models for Ha profile s for substellar objects (10000-12000 K, 
Muzcrolle et al . 2003). This object also has an upper limit o n 
the accretion rate of < 5- 10" 11 M Q yr" 1 (|Gatti et alj|2008h . 
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illustrating that large-scale variability due to hotspots can 
be generated even at very low accretion rates. 

For the third highly variable candidate, an object close 
to the substellar boundary, the dominant cause of variabil- 
ity is probably a change in the circumstellar extinction by 
Ay ~ 2 mag. This may be related to inhomogenities in the 
inner disk at radii 0.01-0.04 AU, which is close to the inner 
edge of the disk. SED modeling confirms an enhanced disk 
scaleheight for this object, which favours the detection of 
disk inhomogenities in lightcurves (Sect.[6]|. 

It is well-established that a substantial number of ac- 
creting objects exhibits h igh-level variability in optical and 
near-infrared bands (e.g. 



1994; Bouvier et al 



.,. Herbst et al . 
19951 ; ICarpenter et al.l l200ll ; lAlves de Oliveira fc Casali 
20081 ). Our results demonstate that this behaviour extends 



into the substellar regime. The fraction of accreting ob- 
jects in the a Ori cluster is 25% for K3-M5 stars and 
14% for objects later than M5, based on the Ha emis- 
sion tlBarrado v Navascues fc Martini [2003). Thus, we ex- 
pect that the total sample of 135 objects published in SE04 
contains about 15-20 accretors (accounting for 30% contam- 
ination) . We find four sources with high-level variability con- 
firmed in at least two independent monitoring campaigns 
(#2, #19, #33, as discussed in this paper, and #43, a likely 
brown dwarf with strong variations in two optical campaigns 
reported in SE04). Thus, the fraction of accreting low-mass 
stars and brown dwarfs exhibiting high-level photometric 
variability is 20-25%. A similar result is obtained when we 
compare the number of highly variable objects in the SE04 
sample (4) with the total number of disk-bearing objects in 
this sample (20, see Appendix |A")) . 

This fact has to be of general concern for studies of 
young stellar and substellar objects: Single-epoch data are of 
limited use for deriving properties of accretors. Probing for 
the presence and extent of variability is important for a reli- 
able assessment of the stellar/circumstellar physics in the T 
Tauri stage. It is usually assumed that accretion affects pre- 
dominantly the wavelength regime < 8000 A through con- 
tinuum excess (veiling) and emission lines. Additionally, 
the disk is expected to produce excess flux at wavelengths 
> 2 /jm. The spectral range around 1 fim, however, is mostly 
assumed to be least affected by those effects and thus rep- 



resents the photospheric fluxes. Therefore, I- and J-band 
fluxes are routinely used to derive s tellar /substellar lumi- 
nosit i es in star forming regions (e.g. iKenvon fc Hartmannl 
1 19951 ; lHillenbrandlll997l ; lLuhman|[l999h . This is also the pre- 
ferred wavelength range for spectral typing and measure- 
ments of effective temperatures. As shown for example in 
this paper, accretion and the presence of a disk can change 
the fluxes in I- and J-band by more than a factor of 2 and 
the J — K colours by more than 20%, thus affecting both 
the total luminosity and the spectral slope around 1 fim in 
a significant way. Variability clearly adds significant uncer- 
tainty to luminosities and effective temperatures for young 
sources, although it is probably not sufficient to explain 
the total amount of scat ter in colour-magnitude diagrams 
l|Burningham et al.ll2005h . 

In this context it is relevant to point out that the masses 
for the targets analysed in this paper have been derived 
in SE04 based on comparison of single-epoch photometry 
with theoretical isochrones. With our new information on 
the origin of the variability, we can improve these esti- 
mates. Since the variability of object #2 is caused by ex- 
tinction, the best approximation for the photospheric flux 
is the lightcurve maximum (minimum extinction). For #19 
and #33 the dominant source of variability is hot spots, thus 
the photospheric flux is best estimated from the lightcurve 
minimum. This gives J-band magnitudes of 14.66, 13.46, 
and 15.33 mag for #2, #19, #33, respectively. Assuming 
a distance o f 420 p c and comparing with the tracks by 
iBaraffe et al.l (| 19981 ) yields approximate masses of 0.09, 0.25, 
and 0.05 Mq. For comparison: SE04 give masses of 0.07, 
0.65, and 0.1 Mq for these three targets. Taking into ac- 
count the extent and the probable cause of variability can 
thus significantly help to constrain the object parameters 
and their uncertainties. 



8 SUMMARY 

In this paper we present a comprehensive analysis of the 
properties of disks around low-mass stars and brown dwarfs 
in the young open cluster a Ori. In the following, we will 
summarize our main results. 
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(i) The typical variability characteristics established of 
solar-mass T Tauri stars are also observed in very low mass 
stars and brown dwarfs. 

(ii) We identify three objects with large-scale variability 
(J-band amplitudes 0.5 mag) sustained over timescales of 
several years, including two possible brown dwarfs. In two 
cases (#19 and #33), the dominant cause of the variations 
are hot spots co-rotating with the objects, in the third case 
(#2) it is variable extinction, caused by inhomogenities in 
the inner disk. The hot spot temperatures are found to be 
significantly lower than 10000 K. The available long-term 
data for these objects make them excellent test cases for 
models of inner disk evolution and magnetospheric accre- 
tion. 

(iii) For accreting objects we caution against using single- 
epoch data in the optical/near-infrared to determine lumi- 
nosities and effective temperatures. At least in some cases, 
the brightness in I- and J-band can change by more than 
0.5 mag, even at low accretion rates. Variability information 
can significantly improve estimates of fundamental parame- 
ters. 

(iv) Combining our dataset with the published results for 
previous seasons, we identify persistent periodic components 
in the lightcurve, covering timescales from ~ 0.5 to 8 days. 
These periodicities may be related to the rotation periods 
of the targets. 

(v) Based on Spitzer data, we find that the three highly 
variable objects all show significant mid-infrared excess con- 
sistent with emission from a disk. Object #2 shows evidence 
for enhanced flaring, which favours the detection of inho- 
mogenities in lightcurves, consistent with the variability in- 
terpretation. 

(vi) We identify 20 objects in a Ori with mid-infrared ex- 
cess due to the presence of a dusty disk, among them 7 with 
likely masses below 0.1 Mq. The typical mid-infrared SED of 
disk-bearing objects is not a strong function of object mass 
for the regime 0.05-0.7 Mq. Thus, evolutionary processes in 
the disks occur do not strongly depend on the mass of the 
central object. 

(vii) Comparing the typical SEDs of very low mass 
sources at 2, 3 and 5Myr, we see a clear trend of declin- 
ing mid-infrared flux levels, best explained by a progression 
in the degree of dust settling. This evolution occurs inside 
out, i.e. it affects the inner disk first and then progresses to 
larger radii. The disk evolution at very low and substellar 
masses is compatible with the current consensus view for 
stars. 
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APPENDIX A: DISKS IN THE SE04 SAMPLE 

The a Ori cluster has been surveyed by Spitzer as part of the 
GO programs #37 (IRAC, PI G. Fazio) and #58 (MIPS, 



PI G. Rieke). These im ages have already been analysed by 
IHernandez et al.l l|2007i 'l. Their work, however, does not fully 
make use of this dataset; for example, it does not give the 
photometry for a significant fraction of the SE04 candidates 
covered by Spitzer. In the following we therefore provide a 
report about the Spitzer photometry for the SE04 candi- 
dates, intended to be complementary to the characterisa- 
tion o f cluster members done by SE04 and IHernandez et al.l 
(2007). The mid-infrared fluxes for the variable sources #2, 
#19, #33 have been used in Sect.[6]for a detailed SED anal- 
ysis. 

Al Spitzer photometry 

We used the pipeline-reduced PBCD images from AORs 
3959552 and 4322304. The IRAC fields cover 94 of the 135 
objects in SE04; 60 of them in all four IRAC channels. The 
MIPS field contains 63 of the SE04 objects (almost all of 
them covered by IRAC, too). The table of photometrically 
selected a Ori members provided bv lHernandez et al.l (|2007f ) 
has only 45 objects in common with SE04. Further 4 of the 
SE04 candidates are contained in their sample of 'uncer- 
tain members', suspected to have lower probability of cluster 
membership. 

Our photometry is based on daophot routines within 
IRAF. Aperture photometry was carried out with a constant 
aperture radius for all objects. For the IRAC images, we 
used a radius of 5pix and a sky annulus of 5-10 pix; for 
MIPS the radius was 5 pix with a sky annulus of 8-13 pix. 
Aperture corrections and zeropoints were applied as given in 
the Spitzer data handbooks. The photometric uncertainties 
were derived by combining Poissonian errors for the source 
fluxes, the error in the sky background estimate, as well as 
~ 5% uncertainty in calibration and aperture correction. 

We compa red our Spitzer photome try with the magni- 
tudes listed in IHernandez et all ((2007) for the 45 objects 
in common. Both datasets agree fairly well. In all four 
IRAC channels; the la deviation in the magnitudes is 0.06- 
0.08 mag with maximum deviations of 0.1 mag (channels 1 
and 2) and 0.2 mag (channels 3 and 4). These differences 
are fully consistent with the uncertainties in the photome- 
try. For the 15 objects with published MIPS fluxes at 24 /im 
the deviations are all within 0.2 mag, again within the er- 
rorbars. 

Only very few of the SE04 candidates have been 
followed-up spectroscopically, i.e. the sample is likely con- 
taminated by field stars in the fore- and background of the 
cluster. The fraction of contaminating objects has been es- 
timated to be ~ 30% (SE04). 

Based on colour-magnitude and colour-colour diagrams 
constructed from Spitzer photometry, we identify disk- 
bearing objects in the SE04 sample. In a first step, we plot 
the 'standard' IRAC colour-colour diagram in Fig. I All 47 
objects have photometry in all four IRAC channels, from 
which 19 fall in the colour range expected for CTTS (dashed 
lines, see I Allen etalll2004 ). The two objects #45 and #116 
have the largest excess colours with [5.8]- [8.0] > 1.5. 

To account for the objects only observed in IRAC chan- 
nels 1 and 3, we plot a colour-colour diagram with [3.6]- [5.8] 
vs. I — J, a proxy for effective temperature. As can be seen in 
Fig. IA2I the 76 objects roughly fall in two groups, which we 
identify as CTTS and WTTS (plus contamination). The 18 
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Figure Al. IRAC colour-colour diagram for the 47 SE04 objects 
with photometry in all four IRAC channels. The dashed lines 
show the limits between CT TS (upper right corner) and WTTS, 
following lAllen et al.l j2004 \. The square marks the object #33. 
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Figure A2. Colour-colour diagram for the 76 SE04 objects with 
photometry in IRAC bands 1 and 3. The dashed line shows the as- 
sumed limit between CTTS (upper part) and WTTS. The squares 
mark the objects #2, #17, #19, #33, and #43. 



objects with [3.6]-[5.8]> 0.5 are good candidates for CTTS. 
In Fig. IA2I we mark objects with squares that have been 
found to be relevant in the monitoring campaigns: #2, #17, 
#19, #33, #43 (the last one is seen as highly variable ob- 
ject in the campaigns by SE04). With the exception of #17, 
these objects are clearly in the CTTS regime, confirming 
that their variability is related to the presence of a disk (see 
Sect. ©. In Fig. QO] the objects #37 and #45 stand out 
with the highest IRAC colours of > 1.0 mag. 

In a third diagram we probe for 24 /im emission for ob- 
jects covered by the MIPS observations. In Fig. IA3I we plot 
K-[24] vs. [3.6]-[5.8] colour for 24 objects with reliable de- 
tection in MIPS channel 1. The horizontal dashed line delin- 
eates the limit between optic ally thick T Tauri disks a nd op- 
tically thin disks, as given bv lHernandez et alj l|2007h based 
on the colours of the /3Pic debris disk. The vertical line 
shows our CTTS/WTTS limit from Fig. QO] 15 objects are 
in the CTTS regime, including the variable ones #2, #19, 
#33, #43. Five more objects show excess at 24 /j,m, but not 
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Figure A3. Colour-colour diagram for the 24 SE04 objects with 
photometry in MIPS band 1 (24 /im). The dashed horizontal line 
shows the assumed limit between optically thick (upper part) 
and optically thin disks, following [Hernandez et al. (2007). The 
vertical dashed line is the limit between CTTS (right) and WTTS, 
as in Fig. EH The squares mark the objects #2, #17, #19, #33, 
and #43. 



at IRAC wavelengths (#17, #21, #26, #93, #109). The 
objects #45, #37, and #116 exhibit the largest values in 
K-[24]. 

For all disk candidates, we re-evaluated the multi- 
wavelength imaging from Spitzer and in questionable cases 
also from SE04 and 2MASS, to exclude possible contam- 
ination from closely neighboured objects. Most objects are 
well-isolated in all images. Objects #26 and #45, both close 
to the detection limit of 2MASS, are likely to be misidenti- 
fications: Instead of the faint source seen in I-band images, 
the 2MASS and Spitzer photometry is dominated by the flux 
from a brighter neighbour at 2-3'.' distance. 

Object #116 has a visual companion at 4" distance 
which contaminates the I-band flux. This star, however, is 
bluer than our target and therefore becomes insignificant 
towards longer wavelengths. The photometry published in 
SE04 probably overestimates the I-band flux and thus un- 
derestimates the I — J colour. Taking this into account, the 
SED still fits the criterium fo r cluster membership (an d was 
also classified as member by Hernand ez et alj ([2007)), but 
I- and J-band fluxes should be treated with caution. 

With Fig. IA1IIA3I we have carried out three tests for 
T Tauri like mid-infrared excess, indicating the presence of 
a disk. The total number of sources that passes at least 
two out of three tests is 22, from which we exclude two 
based on by-eye inspection of the images (see above). Four 
sources from this s ample are investigated in a recent pa- 
per by iGatti et al.l (|2008l ). confirming their status as disk- 
bearing objects. The disk fraction in the sample of SE04 
objects covered by Spitzer observations is thus 20/94 or 
21 ± 4%. 

This can be used to obtain an independent estimate of 
the contamination rate in the SE04 sample. We assume that 
objects with Spitzer excess are confirmed as young cluster 
members, while objects without disk excess may still con- 
tain contaminating objects in the fore- or background of the 
cluster. If we further assume that the actual disk fraction 
among low-mass a Ori members is 34 ± 3%, as derived by 



14 



Scholz et al. 







-1.5 



- A^tJi 






ill 


- 






a 


□ 


- 






i 





□ I 




A 


* 

A 


A 


A ~ 








ill 





1 

Wavelength (/^m) 



Figure A4. Spectral energy distributions for objects with disks 
in the a Ori cluster. All fluxes are normalized to the J-band flux. 
Three objects with unusual SEDs are marked (triangles for #21, 
squares for #37, hexagons for #116). The solid line shows the 
median SED for the high-mass bin, the dashed line for the low- 
mass bin (M <0.2M Q ). 

IHernandez et all l|2007l ). we would expect that the part of 
the SE04 sample covered by Spitzer (94 objects) contains 
59 cluster members (20 x 0.34 -1 ). This would put the con- 
tamination rate at ~ 37%, in agreement with the previous 
estimate in SE04. 

A2 Disk SEDs as a function of mass and age 

The spectral energy distributions (SED) in the wavelength 
range 3-24 [im are mostly determined by the geometry and 
structure of the dust in the inner disk. The most relevant 
parameters affecting the shape of the mid-infrared SED are 
disk flaring and inclination as well as the presence and size of 
inner holes. Dust settling to the disk midplane, thought to be 
accompagnied by grain growth, and inner disk clearing can 
be pro bed by analysing this part of the SED l|Wood et al.l 
12002 j ). 

In Sect. |A"T1 we selected 20 objects in the SE04 sample 
with infrared colours consistent with the presence of a disk. 
18 of them are covered by MIPS, thus have a 24 [im data- 
point. In Fig. IA4I we show the SEDs for these 18 sources; 
all fluxes are plotted on a logarithmic scale and are normal- 
ized to the J-band flux, to facilitate a comparison for objects 
with different photospheric fluxes. 

To look for mass effects in the disk properties, we di- 
vided the sample of 18 disks in two equally sized bins, 
J > 13.5 mag and J < 13.5 mag. This corresponds to a 
separation in two mass bins at roughly 0.2 Mq. In Fig. IA4I 
we overplot the median SED for the high-mass bin (solid 
line) and the low-mass bin (dashed line). As can be seen in 
the plot, the typical SED does not change significantly with 
mass in the range 0.05-0.7 Mq, indicating that the evolu- 
tion of the inner disk in substellar objects occurs on the 
same timescales as in stars. A similar result has been ob- 
tain ed recently for olde r objects in Upper Scorpius (see Fig. 
1 in lScholz et alll2007l ). 

The datapoints for three unusual SEDs are marked with 
special symbols in Fig. IA4I Objects #37 and #116 deter- 
mine the upper limits at all mid-infrared wavelengths. #21 



and #116 show a clear break in the SED around 8 /im, with 
increasing excess towards longer wavelengths, which can be 
interpreted as a possible opacity hole in the disk. 

We probe the evolution of brown dwarf disks by com- 
paring the SEDs of VLM objects in a Ori with younger and 
older objects. For this purpose, we use from our disk sam- 
ple in a Ori the 7 objects with J > 14 mag, correspond- 
ing to M < O.I Mq. To enlarge t he sam ple size, we add 
13 objects from IHernandez et all (|2007l ) with 24 /im de- 
tection and J > 14 mag. As comparison samples, we use 
brown dwarfs w ith 24 fim detectio n in Upper Scorpius (Up- 
Sco. 13 objects. [Scholz et al.ll2007l) a nd Taurus (11 objects, 
IScholz et ai1l2006l ; lGuieu et al.ll2007l ). In terms of their ages, 
these two regions bracket a Ori, with approximate ages of 
2Myr for Taurus and 5Myr for UpSco. 

In Fig. IA5I we show the median SED for all three re- 
gions. In the older UpSco region, the mid-infrared flux lev- 
els are lower than in Taurus, by 0.5-1 order of magnitude. 
The SED in a Ori follows the UpSco SED until 8 fim, but 
is intermediate between Taurus and UpSco at 24 /im. The 
best explanations for lower flux levels in the mid-infrared 
are the onset of inner disk clearing or a lower degree of flar- 
ing, or conversely a higher degree of dust settling to the disk 
midplane. Thus, Fig. IA5I shows a progression of inner disk 
evolution with age. Additionally, the SED in a Ori indicates 
that the onset of evolutionary effects is a function of wave- 
length; it occurs earlier at 8 fim compared with 24 fim. This 
implies that disk evolution occurs inside out, i.e. it begins 
close to the central object and progresses to larger radii. 

A second age effect can be seen by looking at the scat- 
ter instead of the median SED: The total spread/rms of 
the 24 fim fluxes is 3.2/0.8 in Taurus, 1.3/0.3 in <rOri, and 
0.7/0.2 in UpSco. A similar effect can be seen at 5.8 and 
8 fim when comparing Taurus and a Ori (we do not have 
IRAC data for the UpSco objects). Thus, with increasing 
age the scatter in the SEDs decreases, the SEDs and thus 
the inner disks become more similar. 

In summary, the analysis confirms that the standard 
evolution ary picture for T Ta u ri dis ks, as outlined for ex- 
ample in ISicilia-Aguilar et al.l l|2006h . also applies to very 
low mass stars and brown dwarfs. Very young objects show a 
high degree of diversity in their inner disk properties, includ- 
ing a large fraction of highly flared disks. With progressing 
age, dust settling and inner disk clearing affect an increasing 
fraction of objects. This process begins in the disk regions 
very close to the star (~ 1 AU) and progresses to larger radii. 
In our sample in a Ori the effects of disk evolution are not 
a strong function of objects mass. 
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